Analytical contradictions of the 'fixed - node' density matrices 
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Over the last decades the 'fixed-node method' has been widely used for a numerical treatment 
of thermodynamic properties of strongly correlated Fermi systems. In this work correctness of the 
'fixed -node method' for ideal Fermi systems has been analytically analyzed. Rigorous consideration 
shows that the 'fixed-node' calculation of the density matrix even for two ideal fermions leads to 
contradictions. Analogous contradictions results from the virial decomposition of the many fermion 
'fixed - node' density. Numerical results of the 'direct path integral Monte Carlo simulations' show 
f— «) I that the 'fixed - node method' describes the thermodynamic properties of the strongly coupled 

^vj , fermions rather well only at weak degeneracy. Difference in results obtained by these methods in- 

' creases systematically with the growth of the degeneracy at high density and low temperatures. The 

f>, ' reason of this difference results from the replacement of the initial condition by the zero boundary 

^^ , conditions for the density matrix in the Bloch equation in the 'fixed-node approach'. This replace- 

^^ ■ ment results in uncontrolled errors in calculations of thermodynamic quantities due to the wrong 

fT^ ' description of statistical effects. The main conclusion of this work is that the 'fixed-node method' 

^SJ , can not reproduce density matrices of fermions and should be considered as uncontrolled empirical 

approach in treatment of thermodynamics of Fermi systems. 
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CLc I. INTRODUCTION 

. Sh ' Over the last decades significant progress has been observed in theoretical studies of thermodynamic properties of 

C^ . strongly correlated systems of fermions at non-zero temperatures, which is mainly conditioned by the application of 

^~»' numerical simulations (see review [l|). The reason for this success is the possibility of an explicit representation of 

the low-temperature density matrix in the form of a finite-dimensional approximation of the Wiener path integrals 

y without any preliminary physical approximations, which require some physical parameters to be small. The main 

difficulty for path integral Monte Carlo studies (PIMC) of Fermi systems results from the requirement of antisym- 

Cn , metrization of the density matrix in the partition function [3| • As a result all thermodynamic quantities are presented 

as the sum of alternating sign terms related to even and odd permutations of particles. So the thermodynamic 

quantities are determined by the small difference of two large numbers, which arc the sums of positive and negative 

^vj ' terms. The numerical calculation in this case is severely hampered. This difficulty is known in the literature as the 

f^ I 'sign problem'. To overcome this issue sonic approaches have been developed, among which the 'fixed-node method' 

[Bj yna] is widely known. 

The main idea of the 'fixed-node method' [l|, |3|-|5| consists in replacement of initial condition by zero boundary 

f^ . conditions for Bloch equation describing fermion density matrices. So instead of initial condition the 'fixed-node' 

density matrix is to take the 'zero conditions' on the boundaries of the domain of the positive values of an additional 

'trial antisymmetric density matrix' presented in the form of path integral. 

The purpose of this work is to give an analytical proof that this replacement results in contradictions in explicit 
calculations of the density matrix even for two ideal fermions. The analogous conclusion have been analytically 
^ ' obtained twelve years ago in [g from virial decomposition of the many fermion 'fixed-node' density matrix'. However 
paper Q is very difficult for understanding as used the Rueele algebraic approach and was missed by the scientific 
community. This is the reason to discuss correctness of the 'fixed - node method' once more using more simple 
mathematical technique. The main result of this work and paper |6| is that the 'fixed-node method' should be 
considered as an uncontrolled empirical approach in treatment of thermodynamics of fermions. 
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II. DENSITY MATRIX BY THE 'FIXED - NODE METHOD' 

Thermodynamic properties of the fermion quantum system at non-zero temperature are defined by the appropriate 
derivatives of the logarithm of the partition function Q^ ~ Tr{p}. Here p ~ exp(— /3iJ) is the density matrix of a 
quantum system of particles with the Hamiltonian H = K + U equal to the sum of kinetic K and potential energy 
U operators, while /3 = l/fc^T . For our purposes it is enough to consider one dimensional (ID) system of two ideal 
fermions, so U = 0, while the kinetic energy operator is the sum of two kinetic energy operators related to each 
particle K = Ki + K2. Density matrix is the solution of the operator Bloch equation 

with the initial condition p|^=o = 1- 

This operator equation in coordinate representation for fermions can be written in the form [3|, |j| 

= -A {R)pf{R, Ro, P) (2) 

with the initial condition 

PFiR,Ro;0) = Y,i-^r'^HR-PRo) (3) 

p 

where R is the set of coordinates of all particles. One of the possible coordinate representation of the fermion density 
matrix looks like 

a 

where the sum is over the complete set of antisymmetric eigenfunctions (f>a{R) oi H. 

Another exact popular coordinate representation of operator p follows from the operator identity e~^^ = 
g-A/SK . . . g-A;3if g-A/SATjQ^ g^^y ^even of oder unity) integer fixed M. Here the r.h.s. contains M iden- 
tical factors with A/3 = P/M . So one can exactly present the ideal density matrix in the form of finite- difference 
expression of the path integral 

Pf{Rm, Ro; /?) = 1^ ^(-1)"^ fdRi--- dRM-1 ( n PiRk-i,Rk] A/3) J p{Rm-uPRm; a/3) (5) 

' p -' \ fc=i / 

where N is the number of fermions. For N = 2 arguments are two dimensional (2D) vectors composed by the 
coordinates of the first and second particle on ID axes AT*^^' and A^^^ and p(_Rfe_i,_Rfe; A/3) are distinguishable 
particle density matrices. Antisymmetry is put in by the antisymmetrization. The sum runs over all permutations 
with parity up acting on indexes of particles. These density matrices are functions of the 2D arguments Rm^Rq on 
(A(^\ A^^^) plane and inverse temperature /3. 

Mathematical meaning of expression ([5]) for density matrix of two particles is illustrated by Fig. [1] where vectors Rk 
are presented by circles (called often as 'beads'), while density matrices p{Rk-i,Rk', A/3) are denoted by segments of 

lines. Sometimes instead of coordinate < A^ , A^ > it is convenient to use coordinates {7^, rjk} defined by expressions 

7fc = 0.5 (xj^'^ +XJ^A and r]k = (x^^^ - X^Y so Rk = ^XJ^\xj^^\ ^ {'yk,Vk}- Modulus of Jacobian related to 

the change of variables of integration in Eq. ^ from the system of coordinates < A^ t^I \ ^'^ ^^^ system {^k, Vk} 
is equal to unity. Action of perturbation P is illustrated in Fig. [1] by the arrow with letter P (see Fig. [ij . For two 
fermions the sum over permutations is reduced to the sum of contributions of identical and non identical permutations 
with opposite signs. The density matrix has the following general properties: 

Pf{R, Ro; /3) = pf{Ro, R; P) = {-IT'' pf{PR, Ro\ P) (6) 

For further comparisons with the 'fixed - node' density matrix let us remind the solution of the Bloch equation ([2]) 
with the initial condition ([3]). Density matrices in ([5]) are well known for ideal particles and can be written in the 
form : 

TT\Rk-Rk-l\^ \ „^^ / 2£|7fc--7fc-i|^\ „„^ ( 7r|?7fc-?7t_i|^ 



exp ^ Ti '-) exp '-^ ^ expl 

p(E,_i,E,.; A/3) = ^ ^ L = ^ ^ L ^ ^-^ ^ (7) 




• ~F ^ x'^' 



FIG. 1: (Color online) Configurational space of two free fermions. Plotted are the paths related to the two possible permutations: 
identical (below line 7) and non-identical (crossing line 7) permutations. 

where A^ = 2Tih?/S.p/m is the thermal wave length related to A/3. The last factor in ^ for permutation P has the 
form 



p(i?A/_i,Pi?A/;A/3) 



exp ( _-JEIlM^^M^ 



cxpi-H^^2i^_2H^)cxp^ -IP^jM-'jM-iP 
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So the exact well known expression for two particle antisymmetrizcd density matrix looks like [3| 



(8) 
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(9) 



with A2 =: 2TTh?l3/m, Pjm = 7m and Pr/M = -?7a/- K {7a/,'7a/} = {7o,?7o} then P7A/ == 70 and Pr]M = -?yo (see 
Fig. p. 

Now let us consider the basic statements of the 'fixed - node method'. The authors of [3, |j| denote the second 
argument Rq of the density matrix as the 'reference point ' and 'the set of points Rt for which there exists a 
continuous 'space-imaginary time' path with pf{R, Rq] ^ > for < t' < /3 the reach of Rq or r(i?o, /3) '. 
For two ideal fermions the reach can be analytically obtained [j, 0] ■ As example the reach (the half plane 77 > 0) for 
?7o > is shown in Fig. [TJ The reach does not depend on temperature. 

The first wrong basic statement of the papers 0, |j| claims that initial condition ^ can be replaced by a zero 
boundary condition on the surface of the reach. Papers 0, |j| claims that 'it follows because the fermion 
density matrix is a unique solution to the Bloch equation ^^ with the zero boundary conditions (see 
Appendix C) '. However the Bloch equation ^ with the zero boundary conditions is linear equation and, besides 
the trivial solution identically equal to zero, has an infinite number of solutions distinguishing by constant factors, 
while the Bloch equation ^ with initial condition ([3]) has really a unique solution. 

Now let us consider general 'fixed - node' solutions of the Bloch equation in the 7 — 77 plane (sec Fig. [ij. The 
authors of [3, 3 claimed that to obtain the 'fixed - node' density matrix 'one simply restricts the paths in 
Eq. d^ to lie in the T{Rq, (3) '. This means that restriction of integration over i?i, • • • , Rm-i in © by the reach 
gives the exact solution of Eq. ^ with initial condition ([3]). The fallaciousness of this statement can be easily proved 
by explicit integration over variables i?i, • • • , Rm-i in the reach for ID two ideal fermions for M = 2 as well as for 
arbitrary integer M. The boundary surface of the 2D reach for two ideal fermions is exactly known and according to 
papers 0, |j| is the line 77 = in Fig.[TJ So on the whole (X^^^X*^^)) plane the 'fixed-node' general solution of Eq. ^ 
looks like (for two ideal fermions) 
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where 0(r]) is the theta function equal to zero for rj < and equal to unity in the opposite case rj > ■ The theta 
functions restrict to the reach the domains of integration. Now assume that all 'basic statements' of papers [3, Q] are 
correct. Then using the mentioned above general properties of the density matrix we have to admit that integration 
over i?2 ■ ■ • , Rm-1 in the reach gives the exact solution of Eq. ([2]) with initial condition ([3]). So the density matrix in 
the 'fixed-node method' pU)) can be transformed to the following integral over the last variable Ri. 

p,[Ru.R.,^) ~-J dR, ^^ V(M-i)A ^''^ 

exp . - ■■'■"" -' 1 - ^^r. . _ .r\PvM-Vl\' 



V 2(M-1 )A'' yl ^^P V 2 
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A 
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Thus instead of the exact density matrix ([9]) we have infinite number of the 'fixed - node' density matrices (|lip . 
These density matrices take the 'zero boundary conditions' on the surface of the reach {i]m = 0) for any M > 2, but 
depend on two arbitrary constants. More over they depend on complementary error functions, (jlip is correct for any 
fixed M >2. 

For further detail analysis of the function ([TT]) let us consider the limit of A/3 — > {M — > oo). Using definition of 
complementary error functions one can transform (jlip to the form: 

cxp(-H£b-^) 
Pf{Rm,Ro':(^) = — 



4A 

(cxpf-i^ 

x[Ci9{mi} 



(cxp i^- 


Tr\riM-'no\^\ 


-cxp(-^^I^i^))(l + s^5n(^o)) 


2A^ J 


(exp [- 


Tr\VM-na\^\ 


A 

-exp(--l^''-''»l^))(l-,s.g».M) 


2A^ J 



c2e(-rM/) — ^ — = — '- — ^ '- ] (12) 



where 



and er/ is the error function. Here sign{r]Q) is equal to plus unity for 770 > and minus unity for r/o < 0. 

One can easily see from P^ that in the limit A/3 — >■ the 'fixed-node' density matrix for two arbitrary constants 
(Ci and C2) takes the 'zero boundary conditions' on the surface of the reach (for tjm = 0). The 'fix- node' density 
matrix P^ differs from exact expressions Eq. © and contrary to it is a non analytical function. More over the 



'fix-node' density matrix p2|) is identically equal to zero if t]q and i]m are lying in opposite half planes (77 > and 
77 < 0) of the 7 — ?? plane, while this is not the case for exact density matrix ^. All these contradictions mean that 
the integration in the reach can not reproduce the exact solution of Eq. ^ with initial condition ([3]) in spite of the 
'second basic statement' of papers [j, |j|- 

An alternative approach for studies Fermi systems without replacement of initial conditions by zero boundary 
conditions for the Bloch equation is known in literature as the direct path integral Monte Carlo simulation (DPIMC) 
(g. IsMllj. In this approach the sum over all permutations is represented identically as a determinant, which can be 
exactly calculated by the direct methods of linear algebra. The accuracy of this approach depends only on the errors 
of the finite-dimensional approximations of the path integrals and can be improved systematically. Comparison with 
results of the DPIMC simulation show that the fixed - node method describes the thermodynamic properties of the 
strongly coupling fermions rather well at weak degeneracy, when the main contribution to the partition function 
comes from the identical permutation [y, lid . Illj . The difference in obtained results increases systematically with 
the growth of the degeneracy at high density and low temperatures [y, QQ, ^^ . The reason of this difference is in 
the replacement of initial condition by the zero boundary conditions for Bloch equation. This replacement results in 
uncontrolled errors in calculations of thermodynamic quantities due to the wrong description of statistical effects in 
the system of degenerate interacting and non interacting fermions . 

Conclusion. 

The 'fixed - node method' can not reproduce correctly even the two fermions density matrix. Analogous conclusion 
for the many particle density matrix of ideal Fermi system have been analytically obtained in [6| from virial decompo- 
sition. So the 'fixed - node method' can not correctly describe the degenerate Fermi systems. The main result of this 
simple work and paper Q is that the 'fixed - node method' should be considered as uncontrolled empirical approach 
in treatment of thermodynamics of degenerate non interacting fermions. Numerical simulations of thermodynamic 
quantities for interacting fermions by the direct path integral Monte Carlo method results in analogous conclusions. 

Acknowledgements 

Author acknowledge stimulating discussions with Profs. P.R. Lcvshov and S.Ya. Bronin. 



[1] 

[2] 
[3] 
[4] 
[5] 
[6] 
[7] 
[8] 

[9] 
[10] 

[11] 



McMahon J M , Morales M A, Pierleoni C, Ceperley D M 2012 Rev. Mod. Phys. 84, 1607 

Feynman R P and Hibbs A R 1965 Quantum Mechanics and Path Integrals (New York: McGraw-Hill). 

Ceperley D 1991 J. Stat. Phys. 63, 1237 

Ceperley D 1992 Phys. Rev. Let. 69, 331 

Militzer B and Pollock R Phys. Rev. E 2000 61, 3470 

Filinov V J. Phys. A: Math. Gen. 2001 34, 1665 

Ruelle D 1969 Statistical Mechanics. Rigorous results (New York: Benjamin-Cummings) 

Egger R, Hausler W, Mak C H and Grabert H 1999 Phys. Rev. Lett. 82, 3320 and references therein 

Imada R J. Phys. Soc. Japan. 1984 53 2861 

Filinov V S, Bonitz M and Fortov V E 2000 JETP Lett. 72, 245 

Filinov V S, Fortov V E, Bonitz M and Kremp D 2000 Phys. Lett. A 274 228 



